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^ ■ Abstract 
O 

A problem of the crossover from percolation to diffusion transport is consid- 
er 

ered. A general scaling theory is proposed. It introduces phenomenologically 
four critical exponents which are connected by two equations. One exponent 
is completely new. It describes the increase of the diffusion below percola- 
te*) ' tion threshold. As an example, an exact solution of one dimensional lattice 

o ' 

I problem is given. In this case the new exponent q = 2. 

Percolation theory is often used to describe transport properties of disordered systems 
| with a large disorder. The typical problems are random mixture of conducting and non- 

conducting elements, or hopping conduction [|TJ. The regular percolation theory assumes 
that the random elements do not change their positions with time so that the percolation 
paths do not change their spatial locations. The low-temperature electron transport is one 
of examples where this assumption is not valid. Due to electron-electron interaction the 
random potential persistently and substantially changes with time fHH, which may affect 
conductivity near the metal-insulator transition. Such problems appear outside solid state 
physics as well. The class of this problem is known as dynamical percolation. They have been 
studied theoretically using effective medium approximation (See e.g. Ref. || and references 
therein). One dimensional problems of this type have been considered for some models 
without any approximations [El]. 

We concentrate here on the case when the diffusion through conducting media is dras- 
tically faster then the fluctuations of the positions of the conducting and non-conducting 
elements. In this case the transition from diffusion to percolation mechanism has all features 
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of the phase transition and it can be characterized by critical exponents. In this paper we 
introduce a set of these exponents and establish relation between them. For illustration we 
present an exact solution of one dimensional lattice problem. 

Thus, we consider here a random mixture of the conducting and non-conducting elements 
which can change there positions but very slowly. In this situation the resulting diffusion (or 
conductivity) of the particles is non-zero at any small fraction p of the conducting elements. 
The mechanism of this diffusion is as follows. A particle can diffuse only in the conducting 
medium. To move from one conducting element to the other it is waiting until another 
conducting element comes to the element where it resides. At this moment of time a particle 
is able to do a next move. This is the simple diffusion which is characterized by some slow 
waiting time t s . The diffusion inside conducting medium is characterized by much faster 
time Tf. When the fraction of conducting elements p becomes close to the percolation 
threshold p c , but it is less than p c the resulting diffusion increases as l/(p c — p) 9 , where q 
is a novel critical exponent. This happens because the conducting clusters become large. 
But since they are disconnected the particle should wait to jump from one of them to the 
other. Finally, when p — p c is positive and it becomes larger than the width of some critical 
region, the diffusion is described by a regular percolation law D ~ {p — PcY and slow motion 
of the conducting elements is not important. To describe this physics one can use the same 
scaling arguments as for the problem of frozen mixture of elements with large and small 
conductivity ||[7|]. Note that the new problem is not equivalent to the old one, so that 
the exponent q might be different. The scaling hypothesis is valid in the proximity of the 
percolation threshold so that \X\/p c <C 1, where X = p — p c . In this region the scaling 
hypothesis can be written in the form 

Here d = D/Df, h = D s /Df, where Df = a 2 /2r/ is the diffusion coefficient in the conducting 
medium, a is a characteristic length which depends on the model, D w = a 2 /2r s , and ip{Z) is 
some analytical function at all real values of its argument, — oo < Z < oo . We assume that 



h < 1 and that *0(O) = 1. 

Eq. (|I|) contains two independent critical exponents. The meaning of the exponent s is 
that d = h s at p = p c . The exponent m describes the width of the critical region \p— p c \ = h m 
between the percolation and diffusion. 

All the critical exponents which describe d(X) can be related to s and m. At X > and 
X ^> h m the slow changes of the percolation paths are not important. Thus, we have the 
diffusion on the percolation cluster. If the system is larger than the correlation radius of the 
percolation cluster, the diffusion coefficient is related to the conductivity by the Einstein 
relation. Then one has d ~ A*, where the exponent t is a usual percolation exponent which 
describes the conductivity above the percolation threshold. To get this result from Eq. ([[]) 
one should assume that m = s/t. At A < and |A| ^> h m we expect that d ~ h/(\X\) q . 
This gives s + mq = 1. Finally we get two connections between four exponents s,q,t,m, 
namely q = — 1) and m = s/t. 

As an example, we consider below an exactly soluble one dimensional site or bond model 
on a lattice. At the beginning we consider the site model. The sites may be white and 
black with the fractions p and 1 — p respectively. The particle may occupy the white sites 
only and it is able to jump at the nearest site, if this site is white, during the time 77. 
The configuration of the white and black sites slowly changes with time. This change can 
be introduced by many ways. We consider the simplest one. After each time interval r s , 
which is called the renewal time, the configuration of all sites completely changes preserving 
the same p, while the particle remains at the same white site. We assume that r s ^> Tf. 
Since p c = 1 in the one- dimensional case, our solution may only illustrate the increase of the 
diffusion coefficient D ~ (1 — p) q and the width of the critical region |A|, where it deviates 
from this law and tends to the value Df — a 2 /2rj which is reached at p — 1. Since the 
diffusion of a particle within the time interval t s is completely independent of the diffusion 
during other time intervals, one has 

^(t) = (n(ti) + r 2 (t 2 ) + ... + r N (t N )) 2 = rliU) + r 2 2 {t 2 ) + ... + r%(t N ) (2) 
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It follows that 



B=lim r!W = r!M, (3) 



where r 2 (r s ) is averaged over all possible initial position of a particle. 

Now our strategy is as follows. We assume first that the clusters of the white sites are not 
very large so that the particle, which change position during the time Tf crosses them to and 
fro many times during the time r a . In this approximation we find diffusion coefficient D\. 
For the larger clusters we use the continuum approximation, which is valid at 1 — p <C 1 only, 
and find This coefficient gives the correct value at p — 1. These two approximations 
have a large region of overlap if 

^ < 1. (4) 

T s 

Matching Di and D 2 in the overlap region we find the result which is exact if Eq. (Q) is 
fulfilled. 

To calculate D\ we assume that the particle crosses each cluster many times during the 
time t s .Then the average quadratic displacement R 2 (S) is independent of t s and it is given 
by the equation 

^) = ^EX>-ao 2 = ^(s 2 -i). (5) 

D n=lfc=l 

The probability that a particle is within a cluster of S white sites is 

w s = {l-pfp s - 1 S. (6) 
By averaging R 2 (S) over all the clusters one gets 

^) = y ((1 - V? pP^SiS* - 1) = (7) 

Thus, 



where D s = a 2 /2r s . Note that Eq. (|j) has been obtained by Druger et al. ||. The Eq.( 



is valid if r 2 (r s ) << D f r s . It is fulfilled if 



77^ « -• (9) 

At p 1 one gets 

Di = (10) 

r s 

To calculate -D2 we use the continuous approximation which is valid when the size of white 
clusters 5 > 1. This is true if 1 — p <C 1- In this approximation one should solve the 
diffusion equation 

0u(x,t) <9 2 w(x,t) 

assuming zero current at the beginning and at the end of the cluster (x = 0, L) 

du(x, t) 



dx 

The initial condition has a form 



\x=0,L 



0. (12) 



u(x,t)\ t=0 = 5(x - x ), (13) 



where xq is a random point within the interval (0, L). 
The solution has a form 



1 ^ /imx \ / , irnx\ / 7T n ^ \ , . 

w(:M) = - + E cos (-^j cos (—J ex P ( -^r D A ■ ( 14 ) 

The mean squared displacement with respect to the initial position xq is 



L 

2„ 



r2(xo,t)= / (x — Xq) u(x, t)dx 
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(i-^o + ^E^ig + (-ir(l-f)]exp( 



2^2 



7T n 



£ 2 

Finally, taking the average over the initial positions xq on the cluster, we get the time 
dependence of the mean square displacement of the particle on the cluster of the size L. 



At A 3> 1 one can neglect the exponent in the nominator of Eq. (p0|) . Then 

D 2 = D s /(l-p) 2 (21) 

At A <C 1 one can expand this exponent to get D 2 = Df. One can also get the next term 
at small A. Thus, 



Do = Di 



Calculating the integral in Eq. (g2| ) one gets 



1 - ^4 r (~ ~ 1 + exp (-- ) ) dz 

pi A Jo \z V z 



(22) 



D 2 = D f (l - 0.479V 7 !) =D f {l- 1.06^/^(1 - p)) (23) 

at A < 1. 

Note that the linear dependence in 1 — p means that the correction is proportional to the 



ratio of the mean displacement of a particle J DfT s during the time r s to the typical size of 
a white cluster (1 — p)^ 1 . 

The scaling arguments in some different form than Eq. ([]]) can be applied to the one- 
dimensional well. It can be written in a form 

D2 = D f F{^\ (24) 

where h = Tf/r a . It has the same form as Eq. (fj) except that s = 0. This relation is defined 
at p < 1 and 1 -p < 1. In the same way as before, we put F(0) = 1. At large values 
of (1 — p)/h m one has D ~ h/(l — p) q . It follows that mq = 1 which is analogous to the 
relation s + mq = 1 in the many- dimensional case. It follows from Eq.(|2l"D that q = 2. At 
1 — p 1 one can write A = 7r 2 (l — p) 2 r s /2rj. Then it follows that m — 1/2 in agreement 
with scaling relation. The exponent t is not defined in the one-dimensional case. 

Finally we present the solution for the effective diffusion coefficient D which is valid 
at all values of p within the interval < p < 1. It has been shown that D — D\ if 
p/(l— p) 2 Ts/t/- On the other hand D = D 2 if (1 —p) 1. Thus, the two approximation 
have a wide region of overlap (pTf/r s ) - 5 <C (1 - p) < 1. In this region Di = D s p/ (1 — p) 2 



and D2 = D s / (1 — p) 2 . Thus one can get a result which is exact everywhere if r s 3> t/. This 
result is p-D 2 . Finally, 



This result has been derived above for the site problem. One can see, however, that it 
remains unchanged for the bond problem as well. 

The diffusion coefficient D/Df as given by Eq.fl2~5|) is plotted in Fig. 1 versus (1 — p) in 
the double logarithmic scale at different values of T s /rf. The function (1 —p)~ 2 is also shown 
there as a dot-dashed line. One can see that at large values of r s /rf the curves have wide 
regions which are parallel to this line. In these regions the diffusion coefficient increases as 



Finally, we have presented a novel problem which fills a gap between diffusion and per- 
colation in case when the motion of the random media is very slow. We have considered the 
dc transport only. It would be interesting to study a frequency dependent transport in the 
same conditions. 

We thank B. D. Laikhtman for a fruitful discussion. This work was supported by the 
US-Israel Binational Science Foundation, grant 9800097. 
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FIG. 1. The function D/Df is plotted versus 1 — p at different values of t s /tj in the double 
logarithmic scale. The function (1 — p)~ 2 is shown by dot-dashed staight line. 
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